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ABSTRACT 

HD 179949 is an F8V star, orbited by a giant planet at 8 i?^ every 3.092514 days. 
The system was reported to undergo episodes of stellar activity enhancement modu- 
lated by the orbital period, interpreted as caused by Star-Planet Interactions (SPIs). 
One possible cause of SPIs is the large-scale magnetic field of the host star in which 
the close-in giant planet orbits. 

In this paper we present spectropolarimetric observations of HD 179949 during 
two observing campaigns (2009 September and 2007 June). We detect a weak large- 
scale magnetic field of a few Gauss at the surface of the star. The field configuration 
is mainly poloidal at both observing epochs. The star is found to rotate differentially, 
with a surface rotation shear of dO. = 0.216±0.061 rad d~^ , corresponding to equatorial 
and polar rotation periods of 7.62 ± 0.07 and 10.3 ± 0.8 d respectively. The coronal 
field estimated by extrapolating the surface maps resembles a dipole tilted at ^ 70°. 
We also find that the chromospheric activity of HD 179949 is mainly modulated by 
the rotation of the star, with two clear maxima per rotation period as expected from a 
highly tilted magnetosphere. In September 2009, we find that the activity of HD 179949 
shows hints of low amplitude fluctuations with a period close to the beat period of the 
system. 

Key words: stars: magnetic fields - stars: planetary systems - stars: activity - stars: 
individual: HD 179949 - techniques: spectropolarimetry 



1 INTRODUCTION 

Hot Jupiters (HJs) are giant planets, orbiting close to their 
host stars (semi-major axis < 0.1 AU). They represent 
about 25% of all discovered extrasolar planets. Interactions 
between the star and the planet can occur in such sys- 
tems. Star-planet interactions (SPIs) can be of different 
types: tidal interactions due to the proximity and masses 
of the two bodies, and plasma interactions due to the mag- 
netic field, e.g. reconnections between stellar and plane- 
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tary fields (iRubenstein k Schaefe3l2000l: ICuntz et al-llioool : 
iGrieBmeier et al.ll2007l ).l Cuntz et al.l (|20od l suggested that 
such interactions may enhance stellar activity, this enhance- 
ment being modulated by half of the orbital period in the 
case of tidal interactions, and by the orbital p eriod in the 
case o f magnetospheric interactions. Recently, I Fares et"al] 
(I2OIQI ) suggested that the enhancement due to magneto- 
spheric SPIs is more likely to be modulated on the beat 
period (synodic period between the stellar rotation and or- 
bital periods) of the system. 

Observational studies of HJ hosting stars reported that 
not all observed systems show hints of interaction; more- 
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over, for a single system, activity enhancement may be 
present at some epochs, yet absent at other epochs; fi- 
nally, a phase lag between the subplanetary longitude and 
the peak of the activity enha n cement is r eported in some 
systems (jShkolnik et alJ I200I l2005l . I2OO8I ). Different the- 
oretical stud ie s hav e tried to explain these phase lags. 
iMcIvor et alJ (|2006l ) consi dered the cas e whe re the stel- 
lar field is a tilted dipole, IPreusse et alJ (|20Q6h adopted a 
model based on the propagation of Alfven waves within 
the st ellar wind flow relative to the planet and iLanzal 
(|2008l l considered a non- pot ential magnetic f ield configu- 
ration for the closed corona. Cranmer k Saad 11)07) mod- 
elled the Ca II H & K light curve of a HJ hosting star for 
different configurations of the stellar magnetic field (hav- 
ing real solar configurations throughout 11 year cycle), and 
showed that the presence/absence of activity enhancement 
depends on the magnetic field geometry. A statistical survey 
of X-ray emission of stars with HJ suggests that they may 
be ^ 4 times more a ctive than stars with distant planets 
([Kashvap et al. l2008lV but thi s resul t was recently contra- 
dicted by PoDpenhaeger et al.l (|2010l ) who found no signif- 
icant correlations of X-ray luminosity with the planetary 
parameters in their sa mple of stars. Ma^ne tohydrodynamic 
(MHD) simulations of ICohen et"all (l2Q09l ) show that SPI 
may increase the X-ray luminosity. 

HD 179949 is one of the most studied stars for SPI. It 
is an F8 star, orbited by a giant planet with minimum mass 
msini = 0.916 db 0.076M7^, semi-major axis a — 0.0443 d= 
0.0026 AU and orb ital period Porh = 3.092514 ± 0.000032 d 
([Butler et al.ll2006l V The star is reported to undergo epochs 
of activity enhancement modulat ed by the orbital period, 
with a phase lag of - 65° LShkolnik et al.ll2008l ). We or- 
ganised a multi-wavelength campaign to observe this star 
in 2009 September, almost simultaneously in spectroscopy 
(echelle spectrograph @ du Pont telescope @ Las Campanas, 
PHOENIX @ Gemini and SARG @ TNG), spectropolarime- 
try (ESPaDOnS @ CFHT) and X-rays (XMM-Newton). 
With such campaigns one can characterize the system, from 
the stellar surface (activity and magnetic field) to the stel- 
lar corona (X-rays). Simultaneous campaigns also provide an 
accurate description of the large-scale magnetic field allow- 
ing a quantitative modeling of the star-planet interaction. 

In this first paper, we present the results of the spec- 
tropolarimetric campaign of 2009 September and of a previ- 
ous campaign in 2007 June. We describe our data in section 
[21 The magnetic field reconstruction procedure, the mag- 
netic properties and differential rotation of the star are pre- 
sented in section [3l followed by coronal field extrapolation 
and activity analysis in section [H and [S] We summarize the 
results and discuss their implications in section [B] 



2 SPECTROPOLARIMETRIC OBSERVATIONS 

In order to map the surface magnetic field of HD 179949 and 
its evolution, we collected spectropolarimetric data during 
two observing campaigns in 2009 September and 2007 June. 
The data set collected in 2009 September is more exten- 
sive and of better quality than that collected in 2007 June 
- hence we analyze it first in the rest of the study. We used 
ESPaDOnS, a high-resolution spectropolarimeter installed 
at the 3.6-m Canada-France-Hawaii Telescope (CFHT) in 



Hawaii. ESPaDOnS provides spectra that span the whole 
optical domain (370 to 1000 nm) with a resolution of about 
65000 (when used in spectropolarimetric mode). 

The data were reduced using a fully automatic tool, 
called Libre-ESp RIT, installed at t he CFHT for the use 
of the observers (|Donati et al.lll997l ). From collected cali- 
bration exposures and stellar frames, Libre-ESpRIT auto- 
matically extracts wavelength calibrated intensity and po- 
larisation spectra with associated error bars at each wave- 
length pixel. Each spectrum is extracted from four subex- 
posures, taken in different configurations of the polarime- 
ter retarders, in order to perform a full circular polariza- 
tion analysis (|Donati et al.l Il997l ) . Although it is possible 
to extract polarisation spectra from two subexposures, us- 
ing four permits to eliminate (to first order) all systematic 
error or spurious polarisation (as well as producing a null- 
polarisation-test profile, to be discussed below). The spectra 
are normalised to a unit continuum, their wavelength scale 
refers to the Heliocentric rest frame. Telluric lines serve as a 
reference to correct the spectral shifts resulting from instru- 
mental effects (e.g. mechanical flexures, temperature or pres- 
sure variations). The radial velocity (RV) precision of the 
spectra, using this calibration procedure, is about 30 ms~^. 
The RV measurements (listed in Table [1]) are in good agree- 
m ent with the e xpect ations when using the orbital solution 
of iButler et all (I2OO6I ). Fig.[T] represents the RV variations. 

In 2009 September, we collected 19 spectra over 15 
nights (two stellar rotations). Our data have good S/N ratio, 
ranging from 940 to 1910 per 2.6 kms~^ velocity bin around 
700 nm; they sample well the rotation phases. In 2007 June, 
we collected 10 spectra, covering about 1.5 rotational cycles. 
The S/N of the spectra around 700 nm ranges from 950 to 
1280. The complete log of the observations is listed in Table 

m 

The rotational and orbital phases, denoted ERot and 
Eorb, were computed using the two ephemerides: 

To = HJD 2, 451001.51 + 3.092514 Eorb 

To = HJD 2, 451001.51 + 7.6 ERot (1) 

The first ephemeris is that o f iButler et al.l (|2006h , phase zero 
corresponding to the inferior conjunction (i.e. the planet be- 
ing between the star and the observer). For the ephemeris 
giving the rotation phase, we use a rotation period of 7.6 d, 
identified as the equatorial rotation period (see section [3]). 

We apply Least-Squares Deconvolution (LSD) to our 
spect ra in order to improve their S/N ratio (|Donati et al.l 
Il997h . LSD consists of deconvolving each spectrum by a 
line mask, computed using a Kurucz model atmosphere 
with solar abunda nces, temperature of 6250 K (following 
ISantos et al.l |2003 ) and logarithmic gravity of 4.0 cms~^. 
The line mask includes the most moderate to strong lines 
present in the optical domain (those featuring central 
depths larger than 40% of the local continuum, before 
any macroturbulent or rotational broadening, about 4,000 
lines throughout the whole spectral range) but excludes the 
strongest, broadest features, such as Balmer lines, whose 
Zeeman signature is strongly smeared out compared to those 
of narrow lines. The contribution of each line to the average 
line profile is weighted by its depth, so that the strongest 
lines determine the shape of the LSD profile. The cutoff 
value of 40% was selected at the start of the series of stud- 
ies of planet-hosting stars (|Moutou et al.|[2007l : lDonati et all 
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Table 1. Journal of September 2009 and June 2007 observations. Columns 1-8 sequentially list the UT date, the heliocentric Julian date 
(at mid-exposure), the complete exposure time, the peak signal to noise ratio (per 2.6 kms"-*^ velocity bin) of each observation (around 
700 nm), the rotational and orbital cycles (using the ephemeris given by Eq. [T}, the radial velocity (RV) associated with each exposure, 
the rms noise level (relative to the unpolarized continuum level Ic and per 1.8 kms"-*^ velocity bin) in the circular polarization profile 
produced by Least-Squares Deconvolution (LSD) and the false-alarm probability of the detection of the magnetic signature. 



Date (UT) 


HJD 


^exp 


S/N 


Rot. Cycle 


Orb. Cycle 




C^LSD 


fap 


2009 


(2455090+) 


(s) 




(539^) 


(1325+) 


(kms~^) 






25 Sep 


9.73903 


4x1260 


1750 


0.2407 


0.2095 


-24.58 


0.21 


< 10-8 


25 Sep 


9.80652 


4x1190 


1620 


0.2495 


0.2314 


-24.57 


0.24 


< 10-8 


27 Sep 


11.72735 


4x IzoO 


1830 


0.5023 


0.8525 


—24.39 


0.20 


-i n — 06 
10 ^" 


27 Sep 


11.79302 


4x1190 


1780 


0.5109 


0.8737 


-24.39 


0.21 


< 10-8 


28 Sep 


12.73190 


4x1260 


1350 


0.6345 


1.1773 


-24.61 


0.27 


0.08 


28 Sep 


12.79096 


4x1190 


1410 


0.6422 


1.1964 


-24.60 


0.25 


0.008 


29 Sep 


13.73057 


4x1260 


1910 


0.7659 


1.5002 


-24.45 


0.19 


0.0002 


29 Sep 


13.78955 


4x1190 


1870 


0.7736 


1.5193 


-24.44 


0.20 


0.00003 


30 Sep 


14.73560 


4x1260 


1810 


0.8981 


1.8252 


-24.39 


0.21 


0.4 


30 Sep 


14.79547 


4x1190 


1790 


0.9060 


1.8446 


-24.39 


0.21 


0.1 


01 Oct 


15.73477 


4x1260 


1630 


1.0296 


2.1483 


-24.62 


0.23 


0.6 


01 Oct 


15.79425 


4x1190 


1410 


1.0374 


2.1676 


-24.62 


0.28 


0.5 


02 Oct 


16.74999 


4x1260 


940 


1.1632 


2.4766 


-24.56 


0.42 


0.0005 


03 Oct 


17.71701 


2x1260 


980 


1.2904 


2.7893 


-24.40 


0.32 


< 10-^ 


04 Oct 


18.80356 


4x 1260 


1650 


1.4334 


3.1407 


—24.55 


0.23 


< 10-^ 


05 Oct 


19.79055 


4x1190 


1540 


1.5632 


3.4598 


-24.50 


0.25 


0.01 


07 Oct 


21.75237 


4x1260 


1330 


1.8214 


4.0942 


-24.51 


0.27 


0.007 


10 Oct 


24.73615 


4x1260 


1650 


2.2140 


5.0590 


-24.49 


0.23 


< 10-8 


10 Oct 


24.79574 


4x1190 


1480 


2.2218 


5.0783 


-24.51 


0.26 


< 10-8 


2007 


(2454270+) 






(430+) 


(1058+) 








23 June 


4.96537 


4x600 


1280 


0.7178 


0.5095 


-24.49 


0.30 


0.03 


23 June 


5.08479 


4x600 


1260 


0.7335 


0.5481 


-24.44 


0.31 


0.1 


26 June 


8.07869 


4x600 


1060 


1.1275 


1.5162 


-24.45 


0.38 


0.4 


27 June 


8.87788 


4x600 


1200 


1.2326 


1.7746 


-24.36 


0.33 


10-4 


28 June 


9.87806 


4x540 


950 


1.3642 


2.0980 


-24.55 


0.43 


0.9 


01 July 


13.07675 


4x600 


1000 


1.7851 


3.1324 


-24.54 


0.42 


0.2 


02 July 


14.09877 


4x600 


990 


1.9196 


3.4629 


-24.50 


0.43 


0.003 


03 July 


15.00800 


4x600 


1050 


2.0392 


3.7569 


-24.38 


0.38 


0.2 


03 July 


15.07963 


4x600 


900 


2.0486 


3.7800 


-24.37 


0.47 


0.3 


04 July 


15.97775 


4x700 


870 


2.1668 


4.0705 


-24.52 


0.47 


0.5 



l2008l : lFares e t al ." 20091, [ioTol l by optimising SNR in the LSD 
profile as a function of line- dep th cutoff threshold, using the 
method of IShorlin et al In addition to the intensity 

and polarization profiles, LSD produces a null profile (la- 
belled N) used as a polarization check and that should show 
no signal; this helps to confirm that the detected polariza- 
tion is r eal and not due to s purious instrumental or reduction 
effects (jDonati et al.lll997l ). None of the null profiles show a 
polarisation signature. The multiplex gain provided by LSD 
in V and N spectra is of the order of 25 with respect to a 
single line with average magnetic sensitivity, implying noise 
levels as low as 20 parts per million (ppm). 



3 MAGNETIC FIELD AND DIFFERENTIAL 
ROTATION 

3.1 Magnetic modelling 

To reconstruct the magnetic map of the star and estimate 
the surface differential rotation, we use a tomographic imag- 
ing technique, called Zeeman-Doppler Imaging (ZDI). ZDI 
assumes that profile variations are only due to rotational 
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0.0 0.2 0.4 0.6 0.8 1.0 

Orbital Phase 

Figure 1. Radial velocities of HD 179949 derived from Septem- 
ber 2009 (green triangles) and June 2007 (magenta dots) spectra 
as a function of orbital phase, with their error bars. The radial 
velocity curve plotted here i s that using the orbital period and 
amplitude from [Butler et al ] (120061) . 
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modulation and differential rotation, and can invert (in this 
context) series of circular polarization Stokes V profiles into 
the parent magnetic topology and get an estimate of dif- 
ferential rotation at the surface of the star. The magnetic 
field is described by its radial poloidal, non-radial poloidal 
and toroidal components, all expressed in terms of spherical- 
harmonic ex pansions. The expr essions of Br , Be and B(f) can 
be found in lDonati et al.ll2006l : az,m, A,m and 7i,rn are the 
spherical harmonics coefficients describing, respectively, the 
radial poloidal, non-radial poloidal and toroidal components 
of the magnetic field. This description of the magnetic field 
has many advantages: both sim ple and comp lex magnetic 
topologies can be reconstructed (|Donati|[200lh : the energy 
of the axisymmetric and non-axisymmetric modes, as well 
as of the poloidal and toroidal components is calculated di- 
rectly from the coefficients of the spherical harmonics. The 
mean magnetic energy is calculated by the mean over the 
stellar surface of B^ -\- Bq -\- B^. Given the low value of 
the projected rotationa l velocity vsini (7.02 di 0.50 kms~ ^ 
IValenti Fischerll2005l . 6.3 ±0.9 kms'^ lGroot et aLlfTooil ). 
the resolution at the surface of the star is limited. We there- 
fore truncate the spherical-harmonic expansions to modes 
with i < 6. The highest-order spherical harmonics resolv- 
able on a stellar surface typically have orders given by 

i = max (^^^^,4^, where FWHM is the fuh width of 

half maximum of the intrinsic profile. The highest value of 
£ is closely related to the number of surface resolution ele- 
ments around the stellar equator, as described in detail by 
iMorin et all '2010. Sect 3.2. The FWHM of the intrinsic pro- 
file of HD 179949 is 8 kms"^ (see below), while vsmi= 7 
kms~^, implying that modes of order ^ = 5 or 6 are resolv- 
ableEl Whereas the longitudinal resolution depends mainly 
on the vsini, the latitudinal resolution depends also on the 
inclination of the star and the phase coverage of the obser- 
vations. 

To compute synthetic circular polarization profiles, ZDI 
decomposes the surface of the star into 5000 grid cells of 
similar projected areas (at maximum visibility) and calcu- 
lates the contribution of each grid cell to the reconstructed 
profile, given the RV of the cell, the field strength and ori- 
entation, the location of the cell and its projected area. 
Summing the contribution of all grid cells yields the syn- 
thetic profile at a given rotation phase. ZDI proceeds by it- 
eratively comparing the synthetic profiles to the observed 
ones, until they match within the error bars. Since the 
inversion problem is ill-posed, ZDI uses the principles of 
Maxi mum- Entropy image reconstruction (|Skilling BrvanI 
Il984l ) to retrieve the simplest image compatible with the 
data. The form we use for the regularisation function is 
S = - Y^i^^ Ko^lm + Plm + lira)- Morc dctal ls ab o ut th e 
currently use d versio n of ZDI can be found iu iDonatil (|200lh : 
iDonati et al.l (|2006l V Other details and assessment of the 
performance of an earlier version of ZDI code (which was 
not able to recover simple dip o lar magnetic topology) can 
be found in "Brown et al.1 (|l99lh : lDonati Brownl (|l997f ). as 
well as in appendix [B] for the currently used version of ZDI. 

The models we use to describe the local unpolarized 
Stokes I and circular polarized Stokes V profiles at each 

^ We performed a map reconstruction using £max = 4, the results 
are discussed in appendix IA] 



grid cell are quite simple. Stokes I profiles are modelled by 
a Gaussian with FWHM of 8 kms~^and central rest wave- 
length of 550 nm, this Gaussian FWHM fitting best the 
observations. The linear limb-darkening coefficient is set to 
0.6. Stokes V profiles are modelled assuming the weak field 
approximation, i.e. that the local Stokes V profile is pro- 
portional to the effective Lande factor (set to 1.2), the line- 
of-sight projected component of the magnetic field and the 
derivative of the local Stokes I profile. This approximation 
is valid for HD 179949 (whose field strength is a few Gauss). 
The inclination angle o f the star i s approximately 60° given 
the vsini of 7 kms"^ (IValenti Fischer 2 005), the stellar 
radius of 1.19 =b 0.03 Rr;. (jFischer Valentil l2005l ) and the 
rotation period of 7.6 d (see section [3.2. ip . 

Differential rotation (DR) at the surface of the star can 
be estimated using ZDI. Magnetic regions at different lat- 
itudes have different angular velocities, their signatures in 
the spectra repeat with different recurrence rates. We as- 
sume that the rotation at the surface of the star follows 
Q{0) = Qeq — dQsin^{0), where Q{0) and Qeq are respec- 
tively the angular velocities at a latitude and at the equa- 
tor, and dQ is the difference in the rotation rate between the 
pole and the equator. When the data cover more than a ro- 
tational cycle and are well sampled over the rotation cycle, 
one can in principle measure the recurrence rate of magnetic 
regions and thus deduce the differential rotation rate of the 
star. In practice, this is done by reconstructing a magnetic 
map for each pair of differential rotation parameters (Qeq, 
dQ) at a given information content (constant magnetic en- 
ergy), and deriving the associated reduced chi-square x? at 
which modelled spectra fit observations. The optimum DR 
parameters are the ones minimizing Xr- They are obtained 
by fitting the surface of the map with a para boloid around 
the minimum value of x? (jDonati et al.ll2003l V 



3.2 Results 

3.2.1 Differential Rotation 

To determine the DR of HD 179949, we applied the method 
described above to 2009 September data (better sampled 
and covering a longer time span than the 2007 June data 
set). The Xr map we obtain forms a well-defined paraboloid 
shown in Figure [2] The DR parameters we derive are Qeq = 
0.824 ± 0.007 radd-^ and dQ = 0.216 ± 0.061 radd-\ im- 
plying an equatorial rotation period of 7.62 d= 0.07 d and 
a polar rotation period of 10.3 di 0.8 d. The equatorial ro- 
tation period we fin d is compatible with the one found by 
[Shkolnik et al."2003" (< 9 days) but slightly higher than the 
one weakly detected by Wolf & Harmanec 2004 (7.06 days). 
lOjllier CameronI (|2007|) describes dQ as a power- law func- 
tion of the stellar effective temperature; the derived dQ for 
HD 179949 complies with this power law. 

For 2007 June data set, we were not able to measure 
the DR parameters. The data are more noisy than those of 
2009 September and the sampling of the rotational cycle is 
poorer. If we fit the data down to a Xr value of 0.75 (which 
we consider normally as overfitting the data) , we get the DR 
parameters similar to those of 2009 September. 
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n (rad/d) 

Figure 2. Variations of Xr ^ ^ function of fleq, and dfl , derived 
from the modelling of the Stokes V data set for 2009 September. 
The outer colour contour corresponds to a 2.7% increase in the 
and traces a 3 cr interval (99.7% confidence level) for both 
parameters taken as a pair. 



3.2.2 Magnetic configurations 

For 2009 September, the reconstructed Stokes V profiles fit 
well the data down to a x? value of 0.9 (see Fig. [S]). The 
surface magnetic field (Fig. [4j), producing these circular po- 
larization profiles, is mainly poloidal (the poloidal energy 
contributes 90% of the total energy). Low order spherical 
harmonics contribute the most of the poloidal energy, the 
dipolar and quadrupolar contribution being about 40% each. 
40% of the poloidal field is in axisymmetric modes (note that 
non-axisymmetric modes correspond to those with m > 1/2 
and axisymmetric ones to those with m < 1/2). The toroidal 
field constitutes a small fraction of the total field (about 
10%). 

For June 2007, assuming that the DR did not change 
over the two epochs, we fit the observed profiles to a x? 
value of 0.9. The field is weaker than September 2009, but is 
still mainly poloidal (80% of the total energy). The poloidal 
component is mostly axisymmetric (60%). A summary of 
the field properties is presented in Tabled Due to the low 
S/N values we have for June 2007 data, the reconstructed 
magnetic energy could be less than the true energy. A Dis- 
cussion of the effect of S/N and phase sampling on the re- 
construction can be found in lDonati Brownl (|l997l ) and in 
appendix [B] 

Comparing the magnetic field at both epochs, we no- 
tice that it is slightly stronger 1 G) in 2009 September 
than in 2007 June, but this can be due to the better quality 
of the data set. Cool stars having a Rossby number and a 
mass greater than 1 and 0.5 respectively seem to have 
the same trend of magnet ic configurations, i.e. mainl y ax- 
isymmetric poloidal field (jDonati Landstreetl l2009l . Fig- 
ure 3). In this context, the large-scale field configuration of 
HD 179949 is typical for stars with similar masses and rota- 
tion rates. 



4 LARGE-SCALE CORONAL MAGNETIC 
FIELD 

Using the surface magnetic maps, we can extrapolate the 
magnetic field to the corona. In this end, we use the Potential 
Field Source Surface (PFSS) code. T his code was first writ- 
ten t o model the solar coronal field ([van Ballegooiien et al.l 
ll998l V and w as later developed t o be used for stars other 
than the Sun (Ijardine et al.ll2002l ). 

We briefly summarise here the PFSS basic principles 
and underlying assumptions: 

(i) The magnetic fleld B is written as the gradient of a 
scalar potential '0 as B = — V^^, i.e. the field is potential 
(V X B = 0) 

(ii) The three components of the coronal field Br, Be and 
B(f) are described by a spherical harmonics expansion 

(iii) There is a source surface (of radius Rs) beyond which 
the field is purely radial Be{Rs) = Bci){Rs) = 

(iv) The boundary conditions for the radial field at the 
stellar surface are given by the ZDI magnetic maps (when 
having these maps, otherwise the magnetic field is modelled 
by a theoretical model). 

Fig. [5] shows as an illustrative example the coronal field 
lines for the image corresponding to 2009 September, for a 
source surface located at 3.4 R^ (a plausible value given 
that the mean value of source surface radius of the Sun 
is ^ 2.5 Rq). Closed loops, connecting regions of differ- 
ent magnetic polarities, reach different heights in the corona 
(< Rs)- Open field lines have a configuration close to those 
of a dipole tilted by ^ 70° relative to the rotational axis. The 
planet, at 8.1 i?^, is beyond the source surface. Its magnetic 
field can reconnect with open field lines, along which the 
stellar wind can be launched. Given the field geometry and 
the inclination of the star, the open field lines are almost ori- 
ented along the line of sight for some rotational phases (0.2 
to 0.35). The large-scale coronal magnetic field in 2007 June 
is also shown in Figure [S] 

The extrapolation of the coronal field is essential to un- 
derstand the environment in which the HJ evolves. It shows 
for this star that the large scale magnetic field resembles a 
tilted dipole. 



5 STELLAR ACTIVITY 

To study the temporal evolution of chromospheric activity, 
we analyzed two chromospheric activity proxies : Ca ii H & 
K and Ha lines. 

For each observing season (2009 September and 
2007 June) , we calculate a mean profile per proxy. Figure [6] 
shows these mean profiles for 2009 September. The variabil- 
ity in the proxy is calculated by subtracting the mean profile 
to the observed ones. The core of the residual profiles vary 
with time (see Fig.[7|). To estimate the residual emission, we 
fit the line core with a Gaussian (fixing its FWHM and its 
center at the line center); the equivalent width of the Gaus- 
sian gives the residual emission; the residual emission values 
are mentioned in Table (3] The Ca ii H & K and Ha residual 
emission are clearly correlated (Fig. [51 left panel). 

The star exhibits variability on both short and long time 
scales. A modulation of the stellar activity by the orbital pe- 
riod was proposed as a proof of SPI magnetic interactions 
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Figure 3. Circular polarization profiles of HD 179949 for 2009 September (left panel) and 2007 June (right panel). The observed and 
synthetic profiles are shown in black and red respectively. On the left of each profile we show a zbl a error bar, while on the right the 
rotational cycles are indicated. 



Table 2. The field properties for the two epochs of observations. The magnetic field strength averaged over the stellar surface is indicated 
in the second column (the error bar is of ±0.3 Gauss for both seasons). The percentage of the poloidal energy relative to the total one 
and the percentage of the axisymmetric modes (modes with m < 1/2) energy in the poloidal component relative to the poloidal energy 
are reported in column 3 and 4 respectively. The energy in the modes having m = is reported in column 5. 



Epoch 


B 


Epoloidal 


-'--'axisymmetric, m< 1/2 -'--'m=0 




(G) 


% of Etotal 


% of poloidal % of poloidal 


September 2009 


3.7 


90 


36 35 


June 2007 


2.6 


80 


55 54 



(jCuntz et al.ll2000l ). In addition to a modulation with the 
orbital period, a modulation with the synodic period of the 
planet with respect to the stellar rotation is in principle pos- 
sible. The synodic (or beat) period is the period for which 
the planet faces the same configuration of the magnetic field, 
and thus if this configuration is favourable for a magnetic 
interaction (e.g. reconnection) , the interaction will happen 
again after a synodic period. Our aim is thus to search for 
a periodic modulation of the residual activity. We find that 
the activity is modulated, to first order, by the stellar rota- 
tion, and not by the beat period nor by the orbital period, 
see figure [H In order to search for low amplitude periodic 
fluctuation of the activity, two fitting approaches were used 
(a linear and a non- linear one). They will be detailed below. 



5.1 Linear fit to the data 

The open field lines resemble those of a strongly tilted dipole 
(see section m, causing the activity to be mainly modulated 
by half the rotation period. We model the variability by a 
sum of a linear function (for the long term variability) and 
two sine waves (one of period P and one of period P/2). This 
model includes the main frequency of the data and its first 
harmonic (adding a third component at frequency P/3 as in 
iBoisse et al.l 12011 does not modify the result significantly). 
Since the star rotates differentially, the rotational period 
ranges from 7.62 ± 0.07 d to 10.3 ± 0.8 d. For our fit, we fix 
the period and fit the other parameters of the model using 
a rejection criterion to omit the points distant by more than 
10 a to the model. We explore all the domain of rotational 
periods, repeating this procedure for every rotational period; 
the best-fit rotational period is the one minimizing the x? 
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Figure 4. Magnetic maps of HD 179949 for 2009 September (upper panel) and June 2007 (lower panel). The three components of 
the field in spherical coordinates system in a fiattened polar view of the star are presented, down to latitude -30° . The bold circle 
represents the equator. The small radial ticks around the star represent the rotational phases of our observations. The radial, azimuthal 
and meridional fields have the same color scale. 



Figure 5. The extrapolated magnetic field of HD 179949 for 2009 September (left panel) and 2007 June (right panel). White lines 
corresponds to the closed magnetic lines and blue ones to the open field lines (reaching the source surface). The star is shown edged-on. 



of the fit (the number of rejected points never exceeds 2 per 
fit). 

Since we are looking for a possible modulation by the 
beat period and since activity is modulated mainly by rota- 
tion, we first subtract the best fit model to the data and look 
for periodic fluctuations in the residuals (to be called model 
residuals not to be confused with the residual activity). We 
fit these model residuals by a sine wave, and calculate the 
best-fit period of these model residuals. 

For 2009 September, the Ca ii H & K and Ra resid- 
ual emissions are modulated by 7.75lo.55 d and 7.8lJ"^^ d 



respectively, reflecting the fact that both quantities are 
strongly correlated (see figure \T0\ and left panel of figure 
[SI and table [3] for the residual emission values) . The slope 
of the best-fit linear function of Ha vs Ca ii H & K resid- 
ual emission is 0.20. To calculate the false- alarm probability 
(fap) of these signals, we produce 10000 simulated data sets 
by night- shuffling, and fit each data set using the same pro- 
cedure as for our observed data set. The number of data sets 
for which the x? of the best-fit period is inferior to the x? 
of our original signal are divided by 10000 to give the fap 
(Note that we also used the gaussian random noise method 
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Figure 6. Mean Ca ll H &: K (left panel) and Ho; (right panel) profiles for September 2009. 
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Figure 7. The Ca ll H & K residual emission in September 2009 (left panel) and June 2007 (right panel). They show variable features 
at the center of the line. The rotational phases are indicated on the right of each profile. 



to simulate the data sets; the calculated faps are similar to 
the ones calculated using night- shuffling). The fap of the 
Ca II H & K periodicity is about 6%, that of Ha is higher 
(- 30%). 

We then subtracted the best-fit model for each proxy, 
and calculated the model residuals. The Ca ii H & K model 
residuals are modulated by S.S^Q gg d, the Ha ones by 
5.51q;8 d, with a fap of 20 and 7% respectively. These model 
residuals are correlated, as shown in Fig. [8] (right panel). 
The slope of the Ha model residuals vs Ca il H & K model 
residuals is 0.28, higher than that of the residual emission 



(0.20). This correlation may be interpreted as an activity 
enhancement, seen in both proxies, due to the same physi- 
cal phenomenon. The beat period of the system ranges from 
4.27 to 5.23 days, roughly compatible with the period of the 
observed fluctuations (within the error bars). The activity 
enhancement in the model residuals does not correspond to 
a single orbital phase. This enhancement happens at differ- 
ent orbital cycles (0.86, 2.47 and 4.1), contrary to what was 
already reported for this system. 

For 2007 June, the Ca ii H & K residual emission is 
modulated by 8.6l° f d (fap of 0.06%). We do not find a 
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Figure 8. Left panel: Ho; vs Ca ll H & K residual emission for both epochs of observation (2009 September in red and 2007 June in 
blue). The two quantities are correlated, the Pearson correlation coefficient is 0.8. The line is the best-fit linear function to the data, its 
slope is 0.20. Right panel: Ka and Ca ll H &; K residuals for 2009 September after subtracting the best-fit model to the data at Prot- 
The slope of the best-fit linear function is 0.28. 
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Figure 9. The Ca ll H & K residual emission in September 2009 as a function of the rotational phase (left panel) and orbital phase 
(right panel). 
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Figure 10. Ca ll H & K (left) and Ho; (right) residual emissions fitted by a linear function and a sum of two sine waves having as periods 
P and P/2 for September 2009. The best fit period to the data is 7.75 d and 7.8 d respectively. 
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residual emission 


residual emission 




yKIlL S J 


(^Km s J 


9.73903 


-0.007 ±0.040 


-0.002 ±0.007 


9.80652 


0.018 ±0.044 


-0.015 ±0.008 


11.72735 


-0.041 ±0.047 


0.004 ± 0.009 


11.79302 


-0.206 ±0.037 


-0.020 ± 0.008 


12.73190 


-0.285 ±0.039 


-0.010 ±0.005 


12.79096 


-0.224 ±0.028 


-0.019 ±0.008 


13.73057 


-0.110 ±0.038 


-0.012 ±0.006 


13.78955 


-0.011 ±0.035 


0.001 ±0.007 


14.73560 


0.105 ±0.052 


0.018 ±0.010 


14.79547 


0.177 ±0.040 


0.044 ± 0.009 


15.73477 


0.121 ±0.063 


-0.030 ±0.010 


15.79425 


0.066 ±0.039 


-0.027 ±0.007 


16.74999 


0.000 ±0.050 


-0.014 ±0.013 


17.71701 


0.337 ±0.038 


0.101 ±0.014 


18.80356 


0.211 ±0.042 


0.030 ±0.007 


19.79055 


-0.197 ±0.053 


-0.042 ±0.016 


21.75237 


0.310 ±0.045 


0.103 ±0.008 


24.73615 


-0.150 ±0.046 


-0.056 ± 0.009 


24.79574 


-0.137 ±0.037 


-0.072 ± 0.009 


2454270+ 


4.96537 


0.065 ± 0.041 


0.014 ±0.005 


5.08479 


0.053 ± 0.052 


0.014 ±0.001 


8.07869 


0.197 ±0.045 


0.016 ±0.005 


8.87788 


-0.026 ±0.041 


-0.038 ±0.007 


9.87806 


-0.356 ±0.039 


-0.090 ±0.006 


13.07675 


-0.038 ±0.051 


0.009 ±0.009 


14.09877 


-0.044 ±0.051 


0.013 ±0.006 


15.00800 


0.005 ±0.037 


0.011 ±0.004 


15.07963 


0.026 ± 0.049 


-0.016 ±0.006 


15.97775 


0.111 ±0.045 


0.014 ±0.009 



Table 3. Ca ll H K and Ho; residual emission values for Septem- 
ber 2009 (upper part of the table) and June 2007 (lower part of 
the table). 

periodic modulation in the model residuals of the Ca ii H & 
K emission after subtracting the best-fit model. No hints of 
SPI are found for this epoch. 

5.2 Non-linear fit to the data 

We then used a second period search method, consisting of a 
non- linear fit of the data, using our original model : a linear 
function and three sine waves having as periods Pi , Pi/2 and 
P2. We applied this method for the 2009 September data set. 
The periods we find for the Ca 11 H & K residual activity are 
Pi = 8.4 days and P2 = 4.87 days. These periods are com- 
patible with the ones found using the first method within 
the error-bars. For the Ha residual activity, the periods we 
found are 6.0 and 1.27 days. 

Discrepancy between the results found using the two fit- 
ting methods may be due to the limited number of unevenly 
spaced data points, that was also suggested by the high fap 
probability we got using the first method. The source func- 
tion of the Ha line is dominated by the effects of photoion- 
ization and recombination in solar-like stars and F stars, 
so it is not mainly coupled to the local electron tempera- 
ture in the chromosphere but principally to the UV radi- 
ation field of the star whose main component comes from 



the photosphere. The Ca 11 H & K lines are a more suit- 
able proxy for the non-radiative heating produced by mag- 
netic fields in the chromosphere of an F-type star. A definite 
detection/no-detection of low amplitude periodic fluctuation 
needs a larger data set. 



6 SUMMARY 

We present in this paper the first study of the large-scale 
magnetic field of the Hot Jupiter hosting star HD 179949. 
We reconstructed the large-scale surface magnetic field at 
two observing epochs in 2009 September and 2007 June us- 
ing ZDI. The field is mainly poloidal (the poloidal com- 
ponent contributing 80 — 90% of the total energy), with 
a strength of a few Gauss (up to 10 Gauss in some local 
regions). The field configurations, properties, and polarity 
show no main changes when comparing both observing sea- 
sons. We detected differential rotation at the surface of the 
star, the latitudinal shear between the pole and the equator 
being of dQ — 0.216±0.061 radd~^ and the angular velocity 
at the equator being of Qeq = 0.824 ± 0.007 radd"\ The 
equatorial and polar rotation periods are of 7.62 ±0.07 d and 
10.3 ± 0.8 d respectively. The overall field polarity is similar 
during both seasons, thoug h we can not exclude an evolu- 
tion on short timescale (Fares et al.ll2009l V We will continue 
monitoring the system, observing it at least once a year to 
study the evolution of the large-scale magnetic field of the 
star. A comparison between the evolution of the magnetic 
field of similar spectral type stars, with and without HJs, 
is the way to understand whether the presence of HJ influ- 
ences the magnetic field generation (possibly due to tidal 
interactions) . 

This star is reported to show activity enhancement 
modulated by t h e orb i tal pe r iod fo r some observing epochs 
(IShkolnik et al.1 (I200I l2005l . HqqI)). To study possible en- 
hancement due to magnetospheric SPI, we studied the coro- 
nal large-scale magnetic field and the stellar activity. 

The coronal magnetic field is calculated using extrapo- 
lation techniques applied to the surface magnetic maps we 
derived. The field resembles that of a strongly tilted dipole. 

The stellar activity proxies (Ca 11 H & K and Ha ) show 
residual emission modulated by the rotation of the star to 
the first order. We looked for lower amplitude fluctuations by 
subtracting the rotational modulation to the residual emis- 
sion. The new Ca 11 H & K and Ha model residuals are 
periodic and correlated. The variation of these Ha model 
residuals as a function of those of the Ca II H & K shows a 
trend different from the variation of the Ha residual emis- 
sion as a function of that of the Ca 11 H & K (having a 
higher slope value). The periods that best modulate these 
model residuals roughly match the beat period of the sys- 
tem, suggesting that a SPI magnetospheric interaction may 
be at the origin of such an effect. To confirm a SPI activity 
enhancement (or not), a larger data set is needed. 

Further observations of the large scale stellar magnetic 
field are required to study the magnetic cycle length, and 
possible influences of the planet on the stellar magnetic cy- 
cle. Combining observations at different wavelengths will al- 
low us a detailed study of SPI, since it provides information 
for the system from the stellar surface to the corona. 
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APPENDIX A: RECONSTRUCTED MAP 
USING LOWER DEGREES OF SPHERICAL 
HARMONICS 

We have reconstructed the magnetic map of HD 179949 us- 
ing degrees of spherical harmonics up to an order smaller 
than the one used in the paper, for the purpose of compar- 
ison. Indeed, when USmg tmax — 6, only a small amount of 
energy is attributed to orders of ^ > 4 (about 7% of the 
total energy). As shown in the formula used to derive imax, 
it is always possible to resolve up to at least ^ = 4, we there- 
for tested reconstructing the map with imax = 4 for both 
observing epochs. 

The reconstructed map is shown in Fig. lAll Table lAll 
lists the field properties using imax = 4 and imax = 6 for 
comparison. 

We used the same method as in 13.2.11 to quantify the 
DR parameters for 2009 September when using imax = 4. 
We find DR parameters of Qeq = 0.831 =b 0.008 radd"^ 
and dQ — 0.222 zb 0.037 radd~^, implying an equatorial 
rotation period of 7.56 ± 0.07 d and a polar rotation pe- 
riod of 10.3 dz 0.8 d, compatible with the ones found using 
imax = 6. These results are not surprising, since the recon- 
struction procedure uses Maximum- Entropy image recon- 
struction, and thus the simplest image compatible with the 
data is the one to be reconstructed. In our case, low orders 
are sufficient to describe the field in a reliable way. 
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Figure Al. The magnetic map of HD 179949 for 2009 September (upper panel) and 2007 June (lower panel), reconstructed using 
^max = 4. The reader can compare to Fig. [4] where the maps were reconstructed using a higher order of spherical harmonics {imax = 6). 



Table Al. The field properties for the two epochs of observations for the reconstructed maps with imax = 4 and 6. The magnetic field 
strength averaged over the stellar surface (the error bar is of ±0.3 Gauss for both seasons) and the percentage of the poloidal energy 
relative to the total field energy are listed. 
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September 2009 


3.7 
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90 


June 2007 


2.6 


80 


2.9 


83 



APPENDIX B: EFFECT OF LOW QUALITY 
DATA ON THE FIELD RECONSTRUCTION 

In section [3:2:21 we discuss briefly the effect of low quality 
data on th e map reconstruction. S uch effects were discussed 
elsewhere ([Ponati Brownl IiQQTI V for an older version of 
ZDI code. In this appendix, we study the effect of low qual- 
ity data (low S/N, sparse phase sampling), and show that 
the new version of the code (based on spherical harmon- 
ics decomposition) show similar trends as the old one when 
dealing with low quality data. 

Two major cases will be studied, the first one is the 
effect of rotational phase gaps in observation on the relia- 
bility of the reconstructed map, and the second one is the 
reliability of a map reconstructed using a low S/N data set. 



Bl Rotational phase gaps 

To test the effect of rotational phase gaps, we decide to use 
data from 2009 September data set, selecting spectra with 
good S/N ratio. We reconstruct two maps out of two sets of 
data, one with 9 spectra covering rotational phases between 
0.2 and 0.8, and one with 7 spectra covering only half of the 
rotational phases (see Figure IB1|) . Compared to the map 
reconstructed from 19 spectra (section [3231), both of these 
maps show a good recovering of the magnetic field of the 
observed part of the star. Even when half of the rotational 
cycle is not observed, the reconstructed field on the observed 



surface is very similar to the one for the map reconstructed 
with 19 spectra, but features of the unobserved surface are 
missed. 

B2 Low S/N 

One can also question how reliable is a reconstructed map 
based of spectra with low S/N ratio. In order to test this, 
we divide the 19 spectra data set of 2009 September to two 
completely independent subsets, one of 10 spectra with low 
S/N ratio, and a second one with 9 spectra with the highest 
S/N ratio of the original dataset. Table [B2l lists the prop- 
erties of the field in each reconstructed map, i.e. magnetic 
energy and the fraction of the poloidal field. We reconstruct 
the magnetic maps using both subsets, using imax = 6 and 

^max — 4. 

The magnetic energy reconstructed when using a low 
S/N dataset is, as clear from Table [B2l lower than the en- 
ergy reconstructed when using the whole dataset. However, 
the reconstructed maps (not shown here) recover most of the 
reconstructed field of the original map, showing the robust- 
ness of the reconstr uction. These results ar e indeed similar 
to the one found bv lDonati Brownl (|l997l ). 
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Figure Bl. The magnetic map of HD 179949 for 2009 September using 9 spectra with good S/N covering rotational phases from 0.2 
to 0.8 (upper panel) and the reconstructed map for 2009 September using only 7 spectra with good S/N covering half of the rotation of 
the star. The features on the observed surface are well recovered compared to the map in 13.2.21 







^max — 6 




^max — 4 


Number of spectra 


B 

(G) 


Epoloidal 
% of Etotal 


B 

(G) 


Epoloidal 
% of Etotal 


19 (whole dataset) 


3.7 


90 


3.9 


90 


9 spectra with high S/N ratio 


3.7 


90 


3.9 


90 


10 spectra with low S/N ratio 


3.2 


87 


3.3 


87 



Table Bl. The field characteristics when reconstructed using the whole dataset for 2009 September, a subset of only 9 spectra with the 
higher S/N ratio out of the 19 original spectra, and finally a subset of 10 spectra out of the 19 with the lower S/N ratio. The columns 
lists the mean magnetic energy and the energy in the poloidal component of the field, for the reconstructed maps using imax = 6 and 
^macc = 4 respectively (see text for more details). 



